arXiv: 1507.08009vl [cond-mat.mes-hall] 29Jul2015 


Engineering inter-qubit exchange coupling between donor bound electrons in silicon 

Yu E. Wang*, 1 Archana Tankasala, 1 Lloyd C. L. Hollenberg, 2 
Gerhard Klimeck, 1 Michelle Y. Simmons, 3 and Rajib Rahman* 1 

1 Network for Computational Nanotechnology, Purdue University, West Lafayette, IN 47907, USA 

2 Centre for Quantum Computation and Communication Technology, 

School of Physics, University of Melbourne, VIC 3010, Australia 

3 Centre for Quantum Computation and Communication Technology, 

School of Physics, University of New South Wales, Sydney, NSW 2052, Australia 

(Dated: July 30, 2015) 

We investigate the electrical control of the exchange coupling (J) between donor bound electrons 
in silicon with a detuning gate bias, crucial for the implementation of the two-qubit gate in a silicon 
quantum computer. We find the asymmetric 2P-1P system provides a highly tunable exchange- 
curve with mitigated J-oscillation, in which 5 orders of magnitude change in the exchange energy 
can be achieved using a modest range of electric field for 15 nm qubit separation. Compared to the 
barrier gate control of exchange in the Kane qubit, the detuning gate design reduces the demanding 
constraints of precise donor separation, gate width, density and location, as a range of J spanning 
over a few orders of magnitude can be engineered for various donor separations. We have combined 
a large-scale full band atomistic tight-binding method with a full configuration interaction technique 
to capture the full two-electron spectrum of gated donors, providing state-of-the-art calculations of 
exchange energy in 1P-1P and 2P-1P qubits. 

PACS numbers: 71.55.Cn, 03.67.Lx, 85.35.Gv, 71.70.Ej 


I. INTRODUCTION 

Quantum computation has the potential to solve 
classes of problems that are intractable to classical com¬ 
puting, such as large integer factorization which is im¬ 
portant in cryptography 1 . A quantum computer based 
on silicon can take advantage of the existing device fab¬ 
rication platform of the semiconductor industry and the 
exceptional spin coherence times measured to be in ex¬ 
cess of 35 s in recent experiments^. Donor qubits in sil¬ 
icon are promising candidates for a spin based quantum 
com pute r as they have exceptionally long and T 2 

time^lZlSI an d 0 ff er both electron and nuclear spins for 
encoding quantum information 10 utilizing the already 
existing and commonly used silicon device technology. 

While single qubits in silicon have been demo nstr ated 
with both electronic and nuclear spins of donor^ED, the 
next biggest challenge is to demonstrate two-qubit gates 
based on exchange interaction. Ideally, the exchange cou¬ 
pling J in a two-qubit gate needs to be tuned electrically 
by several orders of magnitude between an ’Off’ and an 
’On’ state within a small and realizable bias range. To 
achieve this, the popular Kane architecture uses a J-gate 
between two phosphorus donors to tune the J-coupling 
and A-gates to control individual qubits (Fig. 1(a) (P. 
However, such an architecture leads to a high gate den¬ 
sity, ultra small gate widths, electrical cross-talk between 
gates, and precise donor positioning relative to gates. 
Moreover, the tunability of the exchange coupling is lim¬ 
ited both by the field range the J-gate can produce and 
the field ionization of the electrons to the surface. Calcu¬ 
lations have also shown that the J-coupling oscillates as 


a function of donor separation due to crystal momentum 
state^m, and is therefore sensitive to atomic scale place¬ 
ment errors. All these issues cause severe constraints in 
the implementation of a two-qubit gate in donors. 



FIG. 1: Control of exchange in donor based qubits separated 
by ~15 nm. (a) Kane’s two-qubit donor device with J-gatJ^. 
(b) A MOS device with top detuning gates, (c) An STM pat¬ 
terned device with in-plane detuning gates. (d),(e) schematic 
potential energy profiles of the two-qubit system without and 
with the detuning electric field E, showing a tilt in potential 
energy along the separation direction in (e). 

In this work, we introduce an alternate design for an 
exchange gate in a two-qubit donor system, which al¬ 
lows flexibility in device fabrication and in tuning the 
exchange coupling. In principle, this new design can 1) 
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eliminate the need for additional J-gates between the 
donors, 2) function with a range of donor separations, 
3) provide a 5 order of magnitude J-tunability within 
a modest E-field range of MV/m and lowered ’Off’ 
state exchange, and 4) mitigate the J-oscillations with 
donor separations. This design can also benefit from im¬ 
proved addressability and longer T) timed®, and can be 
used in conjunction with the two-qubit scheme proposed 
in Rel^. This design can therefore overcome some of 
the experimental obstacles for realizing a two-qubit gate 
in silicon. To perform these state-of-the-art calculations 
of exchange energy in donor qubits, we have combined 
atomistic full band electronic structure calculations, elec¬ 
trostatic device simulations, and a two-electron full con¬ 
figuration interaction (FCI) technique over a lattice of 
1.7 million atoms. The calculations therefore describe 
the detailed two-electron spectrum of donor qubits over 
large E-field ranges accounting for both crystal effects 
and electron-electron exchange and correlation effects. 

By studying double quantum dots (DQD) based on 
donor clusters in silicon 13 ", we find that an asymmetric 
2P-1P system outperforms the symmetric 1P-1P system 
in exchange controllability with det uning gates. Anal¬ 
ogous to exchange-tuning in a DQEpEES we envision a 
(1,1) to (2,0) charge transitiorP® in both 1P-1P and 2P- 
1P qubits as a function of a lateral electric field which 
provides the energy detuning. As the electron bound to 
one donor is pulled to the other by the electric field, the 
exchange coupling can be engineered from a small value 
in the (1,1) state to a large value in the (2,0) state due to 
the large spatial overlap of the wavefunctions in the lat¬ 
ter. Such a detuning field can be applied from either top 
gates (see Fig. 1 (b)) in a metal-oxide-semiconductor 
(MOS) device or from in-plane gates realized by scan¬ 
ning tunneling microscope (STM) based lithography 13 . 
Placed on either side of the donor qubits, the detuning 
gates eliminate the need for a sensitive tunnel barrier 
control by the J-gate. Instead, this design realizes a tilt 
in the potential landscape of the two qubits, as shown in 
Fig. 1(d) and (e), and therefore relaxes the more strin¬ 
gent engineering requirements of donor separations and 
gate widths of the Kane architecture, leading to a reduced 
gate density of the computer. 


II. METHODS 

Previous works have calculated the exchange energy 
between two donor-bound electrons in silicon as a func¬ 
tion of separation R using the effective mass approxima¬ 
tion (EMA) and the Heitler-London (HL) formalis m 11 l 17 ( 
which is valid in the regime of small wavefunction over¬ 
lap. Control of exchange with the Kane J-gate was also 
calculated iijniH from EMA based HL method. How¬ 
ever, this method becomes inaccurate at modest gate bi¬ 
ases when the wavefunction overlap increases. The HL 
method also ignores contributions of doubly excited con¬ 
figurations that increase with reduced separations or in¬ 


creasing fields. More recently, the exchange coupling was 
calculated in a 1P-1P system from EMA based molecular 
orbital and configuration interaction approaches for sepa¬ 
ration distances larger than 7.5 nnPOand modest electric 
fieldP! All these works are based on the Kohn-Luttinger 
form of the donor wavefunctions-^ which provides a very 
specific solution to the two-electron problem, and cannot 
provide a full description of the (1,1) to (2,0) charge tran¬ 
sition in which strong Stark effect causes mixing of the 
lowest states with many excited stated 32 -! The atomistic 
configuration interaction method used here goes beyond 
these approximations. 

In this work, we employ a large scale atomistic tight- 
binding (TB) method to simulate the electronic structure 
of P donors in silicon. The full Hamiltonian of silicon and 
donor atoms is represented by the sp 3 d 5 s* atomistic ba¬ 
sis with nearest neighbor interactions 23 ! This method ex¬ 
presses the single electron wavefunction as a linear com¬ 
bination of the atomic orbitals. For a silicon unit cell, 
the model reproduces the full band structure including 
the correct conduction band and valence band extrema, 
effective masses, and band gapPS. A donor is modeled by 
the Coulomb potential of a positive charge screened by 
the silicon dielectric constant, and an onsite cut-off po¬ 
tential Uq adjusted to reproduce the ground state binding 
energy of 45.6 meV below the conduction band for a P 
donor in sP®. It has been shown in earlier works that this 
model captures the central-cell correction, and produces 
correct wavefunction symmetries and energy splittings 
of the donor stateP®. It also captures the valley-orbit 
interaction^, the Stark shift of the donor orbitalP®, and 
the computed ground state wavefunction also showed ex¬ 
cellent agreement with real and momentum space images 
of the donor obtained by STM experiments- 2 ®. Using this 
model, the single electron molecular orbitals of the two 
donors are obtained by, 

# 1,2 = # 0 + H l + H r + eE • r, (1) 

where H\ 2 represents the single electron Hamiltonian 
for electrons 1 or 2, Hq is the crystal Hamiltonian of 
millions of Si atoms, Hl and Hu are the left and right 
donor potential energy respectively, and the last term is 
the potential energy of a detuning electric field E. The 
eigen problem = ed" is solved by a parallel block 

Lanczos algorithm to obtain a set of molecular orbitals 
of the P-P+ system, {dU, 'F 2 , ■•••, ’Pm} corresponding to 
energies {ei, e 2 , e m } where the index goes from the 
ground state 1 to the excited state m. 

Using the m molecular orbitals, all possible two- 
electron Slater determinants (SD) are constructed to rep¬ 
resent non spin-adapted anti-symmetric electron config¬ 
urations of electrons occupying molecular orbitals i and 
j as SDij = ^[^i(ri)^(f 2 ) - 'P J (fi)'F l (r 2 )], where rj 
and r *2 are the coordinates of electron 1 and 2. We em¬ 
ploy a full configuration interaction (FCI) scheme, which 
means that we include all possible single and double ex¬ 
citations of the system for the two-electron problem. We 








3 


found that about 12 single electron molecular orbitals 
(24 spin states) are needed to achieve convergence, which 
leads to the number of SDs being Cf 4 . The two-electron 
Hamiltoinian is given as, 


H 2e = H 1 +H 2 + 


47re(|fi - r 2 1) [fy - f 2 \ 


( 2 ) 


where the last term denotes the Coulomb interaction be¬ 
tween electrons 1 and 2 with the dielectric function being 
e(| -ri — 7 ^ 2 1) ■ If this term is evaluated in the basis of the 
SDs, the familiar Coulomb (C) and exchange (K) inte¬ 
grals are obtained, as follows, 


6 1 

{SDij \^e{\n - r 2 |) In - f 2 | lSDpq) = Gij ’ pq ~ Kij ’ pq ’ 

(3) 

where i,j,p,q < m. The first and the second terms on 
the right hand side are the Cij tPq and Kij, pq integrals: 

c ij,pq = JJ dridr 2 ^*(ri)^*(f 2 ) l _ ^ ^ p (ri)^(f 2 ), 

( 4 ) 

K ij, pq = JJ dndr^* (fi)^*(f 2 ) ^ ^(ri)^ p (r 2 ). 

(5) 

The solution of Eq. [2] in the basis of SDs yields the total 
energies and wavefunctions of the two-electron system, 
which comprises of spin singlets and triplets. The wave- 
functions appear as linear combinations of SDs with coef¬ 
ficients Cj : j , given as T 2f: = JT . C{jSDjj. where i,j < m 
and i ^ j. 

The FCI technique is an exact way to solve the multi¬ 
electron problem only limited by the number of one elec¬ 
tron basis functions used. The method diagonalizes the 
multi-electron Hamiltonian in the basis of all Slater de¬ 
terminants (SD) constructed from the single electron 
states of the system. Each SD represents a multi-electron 
anti-symmetric wavefunction for a particular arrange¬ 
ment of the electrons among the basis orbitals. In ad¬ 
dition to the ground state, the method also captures the 
excited multi-electron spectra, including exchange and 
correlation energieP® This method has been successful 
in solving the challenging problem of the D~ state (the 
two electron state of a single donor) without any fitting 
parameters and providing a charging energy of 45 me\® 
compared to the 44 meV experimental value 31 . The ex¬ 
change energy in the two-electron case can be obtained 
from the difference between the lowest triplet and singlet 
energies, Et and Es respectively, as A Est = Et — Es- 


III. RESULTS AND DISCUSSIONS 

A. Exchange Energy vs inter-qubit separation 

In Fig. 2(b) we show the computed exchange energy 
between two donor bound electrons as a function of donor 
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FIG. 2: Separation dependence of exchange energy, (a) Two 
studied 2P cases with different donor locations on a (001) 
plane in silicon, (b) The exchange energy of 1P-1P and 2P- 
1P qubits in silicon as a function of inter-qubit separation 
along the [100] and [110] crystallographic axes without any 
applied electric fields. 


separation along two high symmetry crystallographic di¬ 
rections [100] and [110]. We consider two cases: 1) a 
symmetric two-qubit system of IP-IP 9 (black triangles 
and squares), and 2) an asym metr ic two-qubit system of 
2P-1P along the lines of RefP^, in which we consider 
two representative locations of the donors inside the 2P 
cluster, shown in Fig. 2(a) as 2P(A) and 2P(B). The ex¬ 
change energy as function of 2P-1P separation for these 
two spatial configurations of the 2P cluster are labeled as 
2P(A) (blue curves) and 2P(B) (magenta curves) in Fig. 
2(b). In all the cases, the exchange energy decreases ex¬ 
ponentially with increasing separation in both the [100] 
and the [110] directions, attaining a few meV for sepa¬ 
ration ~5 nm and a few fxeV or less for ~15 nm. This 
is expected since the wavefunction overlap between the 
qubits decreases as the separation increases. 

The exchange energy of the 2P-1P qubits are found 
to be slightly larger than the 1P-1P qubits even without 
an electric field. This is somewhat counter-intuitive con¬ 
sidering the fact that a 2P system provides a stronger 
confinement for the electron and should minimize the 
overlap with the IP electron. However, due to the asym¬ 
metry of the confinement potential, the 2P cluster is al¬ 
ready slightly detuned towards the (2,0) state with larger 
wavefunction overlap, resulting in a larger exchange. The 
exchange can be reduced below the IP-IP value in this 
case by applying a negative detuning bias, as shown later. 

The exchange energy of the 1P-1P qubits in the [110] 








4 


direction (black triangles) also exhibits oscillations with 
the donor separation. This has been observed in ear¬ 
lier works based on the effective mass HL method^ES 
and has been attributed to the interference between the 
Bloch components of the wavefunctions of the six con¬ 
duction band valleys of silicon. The location of the atom 
sites in the zincblende lattice dictates that the donor sep¬ 
aration in [100] can be increased by at least an unit cell of 
lattice constant ao- However, in [110] the minimum unit 
of separation increases by ao / y/2, which is a separation 
increment of ao/2 in each [100] and [010] directions. 

It is interesting to note that the J-oscillations are 
strongly mitigated in all the 2P-1P cases along both [100] 
and [110] separation. In fact, we considered the two 2P 
cluster configurations, 2P(A) and 2P(B), to verify that 
this is not a one-off case. In 2P(A), the intra-cluster 
P donors are separated by ao times an integer in both 
[100] and [010], and in 2P(B) by ^ times an odd integer. 
These two cases are thus representatives of the donors 
separated by a whole unit cell and a half unit cell. Rela¬ 
tive to this 2P cluster, if the IP distance is varied along 
[100] or [110], the J-oscillations should appear in one case 
in analogy with the 1P-1P system. However, as shown in 
Fig 2(b), we observe strongly suppressed J-oscillations in 
all the 2P-1P cases possibly due to a randomized inter¬ 
ference pattern between the three donor wavefunctions. 
Hence, the 2P-1P unit is a more fault-tolerant system for 
fabrication and control than the 1P-1P unit. 


B. Exchange energy vs detuning electric-field 

Now that we have explored the range of exchange en¬ 
ergies that can be accessed in donor qubits for various 
donor separations, we investigate the detuning E-field 
control of exchange in the 1P-1P and 2P-1P systems. 
Fig. 3 shows the exchange energies as a function of the 
detuning field for 10 nm and 15 nm qubit separations re¬ 
spectively. We consider an electric field range of -2 to 2 
MV/m, consistent with typical E-fields realized in STM 
patterned donor devices. In general, the exchange curves 
show a transition from a low value to a large value for a 
charge transition from (1,1) to (2,0) for the 2P-1P curves. 
This is because when the two electrons are localized on 
separate qubits, the spatial overlap of the wavefunction 
is small, and so is the J coupling. Since a 2P-1P system is 
slightly detuned toward the 2P cluster at zero detuning, 
a negative detuning field further decreases the wavefunc¬ 
tion overlap so that the exchange energy is decreased. 
A positive field pulls both the electrons gradually to one 
qubit, such that the spatial overlap increases, and so does 
the exchange. In the extreme (2,0) regime when both 
electrons are on the left qubit, the detuning field does 
not affect the overlap anymore, and the exchange energy 
levels off. All these regimes are clearly visible in the 2P- 
1P exchange curve for 15 nm separation, in which the 
exchange splitting varies from 48 neV to 5.5 meV over 
a field range of -2 to 0.9 MV/m, resulting in over five 


orders of magnitude tunability. 



Detuning field (MV/m) 


FIG. 3: Electric field induced exchange in donors. Exchange 
energy as a function of a detuning field along the separation 
axis for a 1P-1P and a 2P-1P system with inter-dot separa¬ 
tions of 10 nm and 15 nm. The 2P(B) configuration is used. 

We note that within the 2 MV/m field range, only 
the 2P-1P cases exhibit significant exchange tunability 
in Fig. 3. The 1P-1P case only shows a change from 
0.5 /jeV to 2.4 fieV (5 times) over 2 MV/rn for a 15 nm 
separation distance. Since the 2P cluster has two core 
nuclear charges with one bound electron, it is easier to 
shuttle an extra electron to this system aided by the net 
attractive potential of the core. In the 1P-1P case, the 
electron-electron repulsive energy is stronger due to the 
charge neutrality of each qubit, and a larger E-field is 
needed to reach the (2,0) regime. The calculations show 
that for E< 2 MV/m, the (2,0) regime is never reached 
in the 1P-1P case if the separation is less than 15 nm. In 
addition, since such a system is symmetric, either a posi¬ 
tive or a negative field will increase the exchange energy. 
It can be seen that qubit pairs with larger separations ex¬ 
hibit larger tunability. This is because the same E-field 
causes a larger potential drop between the qubits if their 
separation is larger, and hence provides a larger detun¬ 
ing energy. We also note that the exchange curves show 
that the transition to the (2,0) regime is smoother if the 
separation distance is less. This is due to the stronger 
molecular hybridization of closely spaced qubit pairs. 

After comparing the 1P-1P and 2P-1P cases in Fig. 3, 
one can see that a 2P-1P system with a 15 nm separation 
provides a promising two-qubit unit of a silicon quantum 
computer, yielding 5 orders of magnitude exchange tun¬ 
ability. In our simulations with uniform electric fields, 
we are unable to go to the high field regime for the 1P- 
1P case, as high fields induce a triangular quantum well 
at the lateral domain boundaries causing electron local¬ 
ization in the surface states 22 . However, larger exchange 
energies can possibly be realized even in the 1P-1P case 
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with large spatially varying E-fields from detuning gates 
as shown in Fig. 1(b) and (c). In Appendix, we show that 
the exchange energy of a 1P-1P system can be tuned by 
a factor of 50 using surface detuning gates. We also com¬ 
pare the exchange tunability of corresponding cases un¬ 
der the Kane architecture and show the detuning scheme 
is superior. These simulations include electrostatic sim¬ 
ulations with the Sentaurus technology computer-aided 
design tool (TCAD)^ coupled to the TB-FCI method. 


C. Two-electron spectrum 


We now study the promising 2P-1P system for a 15nm 
separation in more detail including the two-electron en¬ 
ergy spectrum and the singlet and triplet wavefunctions 
at various detuning fields. Fig. 4(a) demonstrates the 
two-electron total energies of several lowest singlet (black 
curves) and triplet (blue curves) states as a function of 
the detuning field. The lowest solid lines are the (1,1) 
and (2,0) singlet states. They mix at low detuning fields 
and anti-cross at 0.7 MV/m, which marks the symmetry 
point of the (1,1) to (2,0) charge transition. The lowest 
(1,1) triplet state represented by the blue dotted line has 
little dependence of the detuning field until 0.9 MV/m, 
indicating that it is decoupled from other triplet states. 
A sharp transition occurs at 0.9 MV/rn, beyond which 
the energies of singlet and triplet states have a linear de¬ 
pendence on the electric field, indicating the occupation 
regime of a (2,0) charge configuration where the lowest 
states of the two qubits are no longer tunnel-coupled. 

Fig. 4 (b) shows the weight of the lowest (2,0) Slater 
Determinant lAot^o-l-) i n the ground state singlet of the 
2P-1P system as a function of the detuning electric field. 
As can be seen, the weight changes from a small proba¬ 
bility 3xl0~ 6 to ~1, indicating a charge transition from 
(1,1) to (2,0). This (2,0) probability can be measured 
by a charge sensor, as in experiments with DQD based 
singlet-triplet qubitJ^EH Thus the same control scheme 
of ReP^ can be utilized to realize singlet-triplet based 
donor qubits with electrical manipulation of J. In such a 
scheme, the system is prepared in the (2,0) state, pulsed 
into the (1,1) state, allowed to evolve through the J- 
coupling, and pulsed back and measured in the (2,0) state 
- all steps being controlled by the detuning field. 

In Fig. 4(c), we show the two-electron density of the 
lowest singlet and triplet states computed from the FCI 
wavefunctions at three different E-fields. At E = 0, the 
singlet and triplet states look similar as they are both in 
the (1,1) charge configuration. The stronger confinement 
in the 2P cluster on the left is responsible for a smaller 
wavefunction extent. At E = 0.7 MV/m, the (1,1) sin¬ 
glet mixes with the (2,0) singlet and the electron density 
gradually shifts to the 2P cluster. However, due to spin 
blockade, the triplet still remains in the (1,1) configura¬ 
tion, with almost negligible change in the wavefunction. 
At high enough detuning fields of 1 MV/m, both the sin¬ 
glet and the triplet are in the (2,0) regime, as verified by 
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FIG. 4: Electric field dependence of singlet and triplet states, 
(a) The total energy of the lowest few singlet and triplet states 
of a 2P-1P system with a 15 nm inter-qubit separation as a 
function of detuning field, (b) Probability of the (2,0) state 
with both electrons on 2P as a function of detuning, (c) 
Two electron density of the lowest singlet (S) and triplet (To) 
states at three detuning field values, showing a transition from 
(1,1) to (2,0). 


the electron densities being localized in the 2P cluster. 

In summary, we have shown that a detuning gate con¬ 
trol of the J-coupling in donors can relax the stringent 
requirements of gate widths, gate densities, and precision 
donor placement needed to realize a two-qubit gate in sil¬ 
icon. In addition, the use of an asymmetric 2P-1P qubit 
pair can yield a giant tunability of the exchange (5 orders 
of magnitude) over a modest field range (3 MV/m) with 
an even lower ’Off’ state exchange than the correspond¬ 
ing 1P-1P qubits. Combined with long Tj times of donor 
clusters^, improved addressability of 2P-1P qubit and 
operation schemes of Ref^, the proposed design helps 
in the experimental realization of the much sought af¬ 
ter two-qubit gate with donors in silicon. The calcula¬ 
tions of detuning controlled exchange energy also helps 
in highly tunable realizing singlet-triplet electronic qubits 
along the lines of DQD^® 
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APPENDIX 

We demonstrate the exchange tunability of the 1P-1P 
and 2P-1P systems in realistic gated device structures as 
shown in Fig. 1. Comparison is made between the Kane 
J-gate (Fig. 1(a)) and the detuning gate structures (Fig. 
1(b) and (c)). For the former case, A represents the hy- 
perfine gates and J the exchange gates, whose width are 
6 nm. The dark red region under these gates is 3 nm 
thick SiC >2 and the light pink region is undoped silicon. 
The qubit-donors are placed right beneath the middle of 
the two A-gates respectively, separated by 15 nm and 
~10 nm below the Si/SiC >2 interface. 0 V is applied on 
the A-gates and the middle J-gate bias is varied from -1 
V to 1 V to tune the exchange energy. The electrostatic 
potential of the device is solved bv a technology com¬ 
puter aided design tool Sentaurus 32 ' and then input to 
the atomistic TB-FCI calculations. 


(a) Kane architecture (b) Detuning gates 



FIG. S: Comparison of tuning exchange energy with Kane’s 
J-gate and detuning gates, (a) Exchange energy as a function 
of J-gate biases for a 1P-1P and a 2P-1P systems with an 
inter-qubit separation of 15 nm along the [100] direction. The 
device structure is shown as Fig. 1(a). (b) Exchange energy 
as a function of top-gate biases for a 1P-1P and a 2P-1P 
systems with an inter-qubit separation of 15 nm along the 
[100] direction. The device structure is shown as Fig. 1(b). 

Fig. S(a) shows the exchange energy of the 1P-1P 
and 2P-1P systems as a function of the J-gate bias. By 
applying negative biases, the exchange energy first de¬ 
creases between 0 V and -0.6 V because the two electron 


are pushed further apart such that their wavefunction 
overlap diminishes. Then the exchange energy reaches a 
plateau between -0.6 V and -1 V because the electronic 
wavefunctions are already pushed as far apart as possible 
by the J-gate bias, resulting in a constant overlap. By 
applying positive biases, the exchange energy increases 
between 0 V and 0.5 V as the two- electron wavefunctions 
are pulled towards the middle region underneath the J- 
gate, and their overlap increases as a result. By further 
increasing the J-gate bias, for the 1P-1P system, the ex¬ 
change energy first decreases (0.6 V) as one of the donor 
electron is ionized to the interface while the other is still 
at the donor site (depending on whether the 6 nm J-gate 
has a bit of asymmetry with respect to the donors). Then 
the exchange energy drastically increases (0.7 V to 1 V) 
because both electrons are ionized to the interface and 
their wavefunction overlap increases sharply. However, 
in this regime, the interface electrons are not distinguish¬ 
able any more and cannot be used as qubit electrons. 

For the 2P-1P system, the exchange energy decreases 
from 0.6 V to 1 V because at 0.6 V, the IP electron is ion¬ 
ized to the interface but the 2P electron is still bound by 
the 2P donors. The interface electron is prone to inter¬ 
face charges and gate noise, which can significantly un¬ 
dermine the coherence time and the qubit performance. 
So this regime for the 2P-1P system is not favorable, ei¬ 
ther. Thus, the exchange can be tuned only by a factor 
of 30 (-1 V to 0.5 V) under Kane’s architecture for both 
the 1P-1P and the 2P-1P systems placed 10 nm deep into 
the lattice. Since the ionization depends on whether the 
surface states near the J-gates are lower in energy than 
the donor states, the ionization bias can vary with donor 
depth. In donors closer to the surface, the electron can 
remain donor bound at larger fields or gate biased? 2 . 

In comparison, we also investigated the exchange tun¬ 
ability with top and in-plane detuning gates. In Fig. 
1 (b), Gi and G 2 are two top gates with a width of 6 nm. 
The qubit donors are placed beneath the center of Gi 
and G 2 , respectively, separated by 15 nm and ^10 nm 
below the Si/SiC >2 interface. Opposite biases with equal 
magnitudes are applied on Gi and G 2 to detune the sys¬ 
tem. Fig. S(b) displays the exchange energy as a function 
of the detuning gate biases for both the 2P-1P and the 
1P-1P systems. For 2P-1P, the Fig. S(b) demonstrates a 
similar exchange tunability, 5 orders of magnitude, to the 
uniform field detuning scheme as shown in Fig. 3 of the 
main text. The in-plane detuning gate structure (Fig. 
1(c)) provides similar tunability as well (not shown). For 
the 1P-1P system, larger spatially varying electric fields 
applied from surface gates increases the J-tunability by a 
factor of 50 (from 0.5 /reV to 22.1 /xeV) as shown in Fig. 
S(b). Compared with a factor of 5 tunability with the 
uniform field detuning scheme, detuning with top-gates 
shows a more promising tunability of the inter-donor ex¬ 
change energy. The exchange tunability in the IP-IP 
system will therefore be limited by the magnitude of the 
electric field realizable in practice. Better performance 
could possibly be achieved with more rigorous device ge- 
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ometry and control optimization. 

In summary, the J-gate in Kane’s architecture provides 
limited exchange tunability for both 1P-1P and 2P-1P 
systems. Detuning gate structures are better designs that 


provide larger exchange tunability. This work therefore 
provides guidance to two-qubit device designs for future 
experiments. 

Electronic address: wangl613@purdue.edu, rrah- 
man@purdue.edu 
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